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Abstract 

We study harmonically trapped ultracold Bose gases with attractive interparticle interactions 
under external rotation in three spatial dimensions and determine the critical value of the attraction 
strength where the gas collapses as a function of the rotation frequency. To this end we examine the 
stationary state in the corotating frame with a many-body approach as well as within the Gross- 
Pitaevskii theory of systems in traps with different anisotropies. In contrast to recently reported 
results [N. A. Jamaludin, N. G. Parker, and A. M. Martin, Phys. Rev. A 77, 051603(R) (2008)], 
we find that the collapse is not postponed in the presence of rotation. Unlike repulsive gases, the 
properties of the attractive system remain practically unchanged under rotation in isotropic and 
slightly anisotropic traps. 
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Bose-Einstein condensates (BECs) have occupied a central role in the study of atomic 
and quantum physics since their first experimental realization. In particular, attractive 
condensates [1], i.e., systems whose bosons attract each other, are distinguished due to their 
peculiar features. Namely, in a trapped three-dimensional attractive gas whose number 
of particles or the strength of the interparticle interaction exceeds a threshold value, the 
kinetic energy cannot balance the (negative) interaction energy and so the gas implodes and 
collapses [2, 3]. However, in the range below the critical interaction strength (or critical 
particle number), there can exist metastable states, i.e., states that will survive the collapse 
for some finite time. The collapse of attractive BECs and situations where it can be hindered 
have been the subject of much interest already two decades ago, see, for instance, Refs. [4, 5]. 
It is furthermore known that the attractive gas, once prepared in a vortex configuration, 
will be more stable against collapse [3, 5, 6]. More recently, fragmented metastable excited 
states [7] and ground states with definite nonzero angular momentum [8] have also been 
found to postpone the collapse. 

On the other hand, it is known that rotating (stirring) condensates is a way of imprinting 
angular momentum in a gas and nucleating vortices [9, 10]. In repulsive gases rotating with 
a frequency smaller than the trapping frequency there can exist configurations where the 
system is well described by a stationary state with some finite nonzero vorticity. Vortices 
[11], vortex lattices [9] and highly correlated - fractional quantum Hall - states [12], as 
well as giant-vortices [13] have all been experimentally observed in repulsive gases. In 
sharp contrast, the behavior of the attractive system under rotation is quite different [14- 
16]. The question of how rotation would affect the stability and collapse of the attractive 
condensate in harmonic traps has recently been addressed [17]. In Ref. [17] it has been found 
at the Gross-Pitaevskii (GP) mean-field (MF) level that the attractive gas can be stabilized 
against collapse for rotation frequencies smaller than the trap frequency. These findings have 
motivated us to attack the same problem at the many-body (MB) level. We show herein that 
rotating an attractive condensate, confined by a harmonic isotropic or slightly anisotropic 
trap, with a frequency below the trap frequency, does not have an impact on the stability as 
well as on the angular momentum of the ground state. We then analyze the problem on the 
GP (MF) level and find as well that no stabilization of the ground state occurs. We stress at 
this point that it has been previously shown that there is no stabilization of the attractive 
gas in an isotropic anharmonic trap with a slight anharmonicity for rotation frequencies 
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below the trapping one (see [14]). 



The system. We consider an attractive BEC of N atoms of mass m, confined by a 
generally anisotropic trapping potential 

V(r) = l -mu 2 [(1 - e)x 2 + (1 + e)y 2 + ( 2 z 2 ] = V (x, y, z) - eV a (x, y), (1) 

where u, e and ( are real nonnegative parameters that determine the frequencies of the trap 
and its deformation, namely cu x = cuy/1 — e, cu y = cuy/1 + e and cu z = tu(, V = ^^-{x 2 + y 2 + 
( 2 z 2 ) is the axially symmetric part of the potential and V a = ^^{x 2 — y 2 ) the 'rotating' 
anisotropy. Since we are interested in the rotating problem we will work in the corotating 
frame of reference, where the MB Hamiltonian takes on the time-independent appearance: 



N r j- 9 -i N 

fir 



2m 



+ \ o yS(r i -r j ), (2) 



i<j 



where Q is the frequency of the rotation around the z-axis, L z the z-projection of the angular 
momentum operator and Ao measures the interaction strength and takes on negative values 
for attraction. We set hereafter h = m = u = 1 so as to work in dimensionless units. 

The Hamiltonian of Eq. (2) admits exact solutions in the absence of interaction, i.e., 
when Ao = 0. In the case of isotropic system (e = 0, ( = 1) and in the limit of fast 
rotation (Q — > u) the energy levels are organized into what is known as Landau Levels. 
The same holds true for weak interparticle interactions [18-20]. Thus, in the fast rotation 
and weak interaction limit the Lowest Landau Level (LLL) is particularly designated for the 
description of the ground state of the system. The orbitals that comprise the (scaled) LLL 
have the form ip^ LL (r) = N k r k e~ r2,/2a2 Y^(9, 0), k = 0,1,2,..., where Y£ is the spherical 
harmonic with I = m\ = k and N k is the normalization constant. The scaling parameter 
a defines the width of the Gaussian part and will be treated variationally, i.e., so as to 
minimize the total energy. Of course, if Ao = then a — 1. At the resonance, Q r = u>, all 
(infinitely many) orbitals of this set become degenerate in energy. The above orbitals can 
also be expressed in Cartesian coordinates as appropriate linear combinations of the solutions 
(f>i(x, y, z) = (p nx (u x , x)(p ny (uy, y){p nz (u z , z) of the three-dimensional harmonic oscillator, i.e., 
the scaled Hermite-Gauss functions, (p nx (u> x ,x) = ~^^j~H nx (^^^x^j e~ UJxx2 ^ 2a2 , where 
H n (. . . ) denotes the Hermite polynomial of degree n. Namely, for the isotropic case u x = 
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ujy — uu z , we rewrite the orbitals as 

with q = (ipl! LL \(f>i) , n x +n y = k, k = 0, 1, 2, . . . , and % = i(n x , n y ) = ^[n x +3n y + (n x +n y ) 2 } is 
a function that we employ to uniquely map the pair {n x , n y } to the single parameter i. Once 
we depart from the isotropy of the trap the infinite degeneracy now at Q r = uj x < u> y , is not 
lifted [21] and the above LLL states are not solutions of the anisotropic system. Since the 
radial symmetry of the trap is broken the orbitals do not possess exact angular symmetries 
and one cannot express the solutions in terms of pure spherical harmonics anymore. Instead, 
one should resort to the orbitals ipk expressed as a mixture of functions 4>i- The same 
transformation coefficients q, that are defined above for the isotropic case, can also be used 
for generic u x ^ u y ^ u z . This transformation maps the functions from the Hermite-Gauss 
representation to that with nonzero (expectation value of) orbital angular momentum. Of 
course, for e = and ( = 1 (i.e., for isotropic traps) the mapped orbitals give back the 
spherical harmonics. The expectation values of L z for the orbitals ipk, for e small enough, 

are {i>k\L z \ipk} = h — 1 + y + 0{e A ) k, with k — 0, 1, 2, Note that, when e ^ 0, the 

orbital set {ipk} is also not an exact solution of the noninteracting anisotropic Hamiltonian, 
since the linear combination ipk = Yl mixes nondegenerate states. However, in the limit 
of small e, this choice is justified on account of working with single-particle states ipk that 
have nonzero (expectation value of) angular momentum Ik and thus allows for a possible 
coupling to the rotation. 

Many-body approach. We study our system at the MB level, i.e., beyond a MF de- 
scription. To this end we follow the Configuration Interaction (CI) expansion, a general 
variational MB method that allows the system to fragment and takes into consideration 
fluctuations of the states. For details on this method and the construction of the configu- 
ration space the reader is referred to the literature, e.g., [8, 22-24]. The MB wave function 
|\&) of the system is expanded over a set of functions (permanents) , 

|*>=X)C7 i |$ i >, (4) 

each describing a MF state of a condensed or fragmented Bose gas of iV atoms. The per- 
manents are built over a certain set of M single-particle functions (orbitals). In this work 



M = 4 and the set of orbitals comprise the LLL and its anisotropic extension, as de- 
scribed above. The permanents can be written in an occupation-number-representation as 
= \n) = \n , rii, . . . um-i), where it is meant that rii bosons occupy the 4>i orbital, 
satisfying £\ = N. The Hamiltonian of the problem is then represented as a matrix % 
over the permanents and diagonalized. The eigenvalues E^ of H are the energies of 
the states. The eigenvectors {Cj} of % provide us with the wave functions with which one 
can compute various quantities like the natural occupation numbers pi of the ground and 
excited states, with X^P« = N. Note that the natural orbitals and the orbitals described 
and used in the expansion above coincide. This holds for the isotropic case (due to the sym- 
metry of the problem) and has been found (numerically) to be well satisfied for the slightly 
anisotropic case discussed below. From the natural occupations we can calculate the total 
angular momentum of the ground state as L = Y^n=o^Pi- By varying the parameter a (i.e., 
the Gaussian width of the orbitals), we minimize the energies per particle e = E/N as a 
function of the rotation frequency Q, for some fixed value of A = |Ao|(iV — 1) and determine 
the optimal value <r . The analysis of the system that follows is always done for optimal 
states, i.e., at er . The number of particles is hereafter set to N = 12. 

We denote with A c the critical value of the parameter A where the ground state of the 
condensate ceases to exist. This is calculated as the largest value of A where there is a (local) 
minimum in the energy E as a function of a. The absence of such a minimum denotes a 
collapsed state (see also [7, 8, 25, 26]). We are interested in the dependence of A c on the 
rotation frequency Q. In Fig. 1 we plot the critical value A c against Q for the isotropic 
e = 0, C — 1 an d the slightly anisotropic case e — 0.1, ( — 1. 1 The values of f2 range 
from to Q r = u x = y/1 — e. At exactly the resonance frequency Q r , the energy diverges 
and the gas becomes mechanically unstable. We notice no change in the stability of the 
ground state of the isotropic system as the rotation frequency Q increases from up to 
and only a negligible increase in A c (f2) of less than 0.1% for e = 0.1. The value of the 
critical parameter A c = 8.425(9) remains unchanged when e = for the whole allowed region 
of Q, and marginally increases from A c (0) = 8.436(2) to A c (fi r ) = 8.440(7) for e = 0.1. 2 

1 Note that even smaller trap anisotropies are sufficient to nucleate vortices in experimental setups, like 
e = 0.025 for instance, in the rotating repulsive gas of Rcf. [10]. 

2 Here and hereafter, when we write X c (fl r ) it is meant, mathematically, A c (0) in the limit of the resonance 
frequency Q — > fl r . The same is meant for other system's properties at the resonance frequency. 
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These results, obtained at the MB level, obviously contradict the GP results of Ref. [17] (see 
analysis and discussion below). 

Next, to analyze the MB results, we chose A = 3 as a representative value of the interaction 
parameter of an isotropic system (e = 0) with noncollapsed ground state and calculated the 
energy per particle e, the angular momentum per particle L/N and the natural occupations 
Pi, i — 0, . . . , 3 for the ground state. We found that the above quantities remain constant 
for any Q G [0,fi r ). The state remains condensed (po = N), carries no angular momentum 
(L/N = 0) and has energy e = 1.396(8). For the anisotropic case of e — 0.1 we also found 
that the above quantities practically do not change. Namely, p marginally decreases from 
Po(0) ~ 12 to po(^r) — 11.998(9), and the rest of the natural occupations change from 
pi(0) ~ l(T 8 ,p 2 (0) ~ 10- 5 ,p 3 (0) ^ 10- 13 to Pl (Q r ) ~ 10~ 7 ,p 2 (Q r ) = 10~ 3 ,p 3 (Q r ) ~ lO" 12 . 
There is an insignificant decrease in the energy [from e(0) = 1.395(8) to e(Q r ) = 1.395(4)] and 
a corresponding increase in the angular momentum [from L(0)/N ~ to L(Q r )/N = 2-10 -4 ]. 

The fact that the ground state of the isotropic system is found to be fully (i.e., 100%) 
condensed deserves some discussion. The total absence of depletion and fluctuations in this 
case is explained if one considers the MB orbital set used: since each orbital ip^ LL has 
different angular symmetry a coupling between the different modes is forbidden due to the 
symmetry of the problem. Any nonzero occupation of the i = 1,2,3 orbitals would result 
in the change of the total angular momentum of the system. Naturally, such a coupling 
is induced in the system when the anisotropy e is turned on and hence the occupations 
Pi,i = 1,2,3 can be nonzero. Nonetheless, as we have found above, for attractive systems 
in three-dimensional isotropic and slightly anisotropic traps, coupling of the ground zero- 
angular-momentum state to excited-states with nonvanishing angular momentum essentially 
does not occur, even for rotation frequencies as high as the resonance frequency Q r . In other 
words, for rotating attractive BECs none of the (or, for slightly anisotropic traps, 

ipk>o) states becomes the state lowest-in-energy, even for rotation frequencies as high as the 
resonance frequency Q r . 

Do the above findings change for a MB basis set that does allow for ground-state de- 
pletion? The answer is negative. Having used, in place of the LLL, the set consisting of 
the s,p + ,po and p_ orbitals (see in this respect Ref. [8]), we found the ground state of the 
isotropic system slightly depleted (i.e., about 98% condensed for A ~ A c ), but its angular 
momentum zero for all rotation frequencies up to the resonance frequency Q r . Side by side, 
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the depletion and fluctuations of the ground state do not depend on the rotation frequency. 
Importantly, the critical value A c for the collapse does not depend on the rotation frequency 
as well. The same conclusion holds for slightly anisotropic traps (e = 0.1). In summary, we 
have shown by a MB approach that the critical value of the interaction for collapse, A c , of 
rotating three-dimensional attractive BECs does not depend on the frequency of rotation. 

The fact that the ground states of both the isotropic (e = 0) and the slightly anisotropic 
system (e — 0.1) were found at the MB level to be essentially fully condensed for any rotation 
frequency Q smaller than the resonance frequency fl r , means that the GP theory should be 
valid here and reproduce the MB conclusions. 

Analysis within the Gross-Pitaevskii approach. We now want to turn from the 
MB to the GP (MF) description and address the same question, namely how the stability 
of the attractive gas is affected as the system is rotated externally. The GP theory assumes 
that all particles reside in the same single-particle state and hence the wave function for 
the state of the whole system is given by a single permanent ^ GP = Y\f i>Gp{ Y i)- The GP 
orbital tpcp f° r the ground state of the rotating gas should be represented with an ansatz 
that takes into consideration orbitals with nonzero angular momentum, as done in the MB 
treatment. To this end we expand ipcp as a linear combination 

fcp(r, cr ) = ^2 & ^( r ' °")> ( 5 ) 
k 

where the basis ipk is the same as the one used in the MB computations reported above. The 
coefficients bk and parameter a (Gaussian width of the orbitals) are determined variationally 
with the normalization constraint \ bk\ 2 = 1 and the summation running over from k = 
to k — 3. We calculate the expectation value E = (^gp\H\^gp) with the above GP ansatz 
and minimize it with respect to the parameters b iy i = 1,2,3 and a, for different values of 
the interaction parameter A = |A |(iV — 1) and for given values of il G [0, y/1 — e) and the 
(small) trap anisotropy e. The expectation value of the angular momentum operator for 
ipGP is I = j £ijb*bj, i,j = 1,...,4, where the matrix elements £y = (ipi\L z \ijjj) are 
given, in second order approximation, as Ca = = ^1 + • (i — 1), £13 = £31 = 
£24 = £42 = 70§ and the rest of the elements are zero. The total angular momentum is 
L = Nl. 

We calculate the critical value of the interaction A c as a function of the rotation frequency 
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Q, for the cases of e — 0, ( = 1 and e — 0.1, C — 1- As anticipated from the MB analysis, we 
again found no essential change in the stability of the gas, as fl varies from to f2 r = y/1 — e. 
Namely, in the isotropic case, the GP ansatz of Eq. (5) yields the value A c = 8.425(9) which 
coincides with that obtained from the MB analysis, and remains fixed for any Q e [0, fl r ). 
In the case of anisotropic trap (e = 0.1) we found A c (0) = 8.436(3) and a negligible increase 
as fl increases, i.e., A c (fi r ) ~ 1.0005 • A c (0). We then fix the interaction parameter to A = 3 
as before. In the isotropic case, e — 0, the energy e = 1.396(8) and angular momentum 
L/N = remain constant for all f2 G [0, Q r ) and, as above, the values coincide with those of 
the MB ansatz. For e = 0.1, we found e(0) = 1.395(8), e(fi r ) = 1.391(8), L(0)/N ~ and 
L(Q r )/N = 0.045(2). Namely, the energies found in the MB and GP approaches are almost 
identical while the angular momentum computed within the GP theory at the resonance 
frequency Q r is somewhat above the value that the MB theory gives. Nonetheless, both 
values of angular momentum can be considered practically zero. In conclusion, the rotation 
does not increase the stability of the ground state described by the GP ansatz of Eq. (5). 

Last, we re-examine the attractive rotating gas using a different GP ansatz that has 
been previously used in the literature, namely the ansatz of Ref. [17] (see also references 
therein). The authors of Ref. [17] considered a GP ansatz for the ground state of the system, 
which they expressed - depending on the geometry of the confining potential - either as a 
Gaussian-sech single-particle wave function: 

0(r) = [N{2l x l y l z -K)- 1 ] 1/2 e- x2 ^ 2 e- y2 / 2l y 2 sech (j-^j e mxy (6) 

or as a Gaussian: 

0(r) = [N(lJ y l z )-\- 3 / 2 ] 1/2 e -* 2 /** 2 e -y 2 Wy 2 e -* 2 Wl e i™y. ( 7 ) 

The parameters l x ,l y ,l z and a were to be determined variationally. The phase axy put in 
Eqs. (6)-(7) is referred to as the 'quadrupolar flow' term; a nonzero value of a increases the 
energy of the isotropic system in this state. 

In Ref. [17] it is found that, using the ansatz of Eq. (6) or (7), the stability of the gas 
is significantly increased with increasing frequency Q. However, an algebraic error in the 
above work is responsible for this (erroneous) behavior of the energy per particle e as a 
function of the rotation frequency f2. Already in Eq. (4) of Ref. [17] there is a sign error; 
redoing carefully the calculations we convinced ourselves that in the true expression the sign 
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in the last term of the integrand is a plus instead of a minus. This sign error gives rise to 
an extra (negative) term in the GP energy functional which overestimates the dependence 
of the energy on Q and artificially reduces the energy of the system (see the Appendix for 
more details). 

In fact, also with the ansatze of Eqs. (6)-(7) the external rotation does not practically 
affect the stability of the condensate, in the sense that the critical value of the interaction 
parameter A c does not essentially change with f2. We have verified this by calculating the 
energies and critical parameters A c by varying all three parameters l x , l y , l z of Eqs. (6)-(7) for 
( = 0, 1, 5 and e = 0, 0.1, as it is originally done in Ref. [17] (the parameter a is expressed 
as a function of l x and l y and absorbed into the GP energy functional as in [17]). The 
Gaussian-sech ansatz of Eq. (6) is used in the case ( = 0, while the Gaussian ansatz of 
Eq. (7) is used when ( = 1 and 5. The critical A c of the radially symmetric systems (i.e., 
e — 0, £ — 0, 1,5) remains fixed, while in the slightly anisotropic systems (i.e., e = 0.1, 
£ = 0,1, 5) A c does not increase more than 0.2% as f2 increases from to Q r . The computed 
values of A c for different values of £ and e are presented in Table I. We then fix A = 3 and 
calculate the energies and angular momenta of the ground state. For e — 0, irrespective of 
the choice of £, the energy has been found to be independent of the rotation frequency Q. 
For instance, we found e(e — 0, £ = 1) = 1.396(8). For e = 0.1 the energy e(Q r ) was found 
to decrease by 1.5% for £ = 0, by 0.32% for £ = 1, and by 0.04% for C = 5 with respect 
to the corresponding energy e(0) of the nonrotating system. Last, the expectation values of 
the angular momentum of the two ansatze are (almost) exactly zero for e = (e = 0.1), 
regardless of the value of £. 

Finally, we point out that, for the cases examined, the optimal value of parameter a is 
practically zero. Note that the minimization of the GP energy functional with the ansatz 
of Eq. (5) yields a distribution for the coefficients 6, (6 — 1 and fej ~ 0, i — 1,2,3) which 
essentially includes only the first of the LLL Gaussian-shaped orbital. I.e., the ansatze of 
Eqs. (7) and (5) essentially coincide with the respective orbital of isotropic systems. 

Summary and Outlook. We have studied the stability under rotation of attractive 
ultracold Bose gases, confined by an isotropic as well as a slightly anisotropic harmonic trap. 
The problem has been mapped to and calculated in the corotating frame. Both many-body 
and Gross-Pitaevskii approaches revealed that the rotation does not affect the stability of 
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the gas against the collapse. Namely, the maximum value of the interaction strength A c 
where the attractive gas collapses remains essentially unchanged as the rotation frequency 
Q varies within the extreme values and Q r = u x = y/1 — e, where tu x is the frequency of 
the trap in the direction of the weakest confinement. 

We have found here on both the MB and the GP (MF) levels that the ground state of 
the rotating attractive system carries zero (or almost zero in the anisotropic case) angular 
momentum for the whole range of the allowed values of Q. Obviously, no vortex states are 
created. In the MB treatment, this means that no transition between LLL states of different 
angular symmetries has been found for the rotating attractive system. In the GP (MF) 
analysis this means that no symmetry broken states were found to be energetically favorable 
as the rotation frequency Q increases from to Q r . In both MB and GP approaches, the 
energy of the ground state remains practically unchanged and the attractive gas is condensed 
in the nodeless s-orbital as the rotation frequency increases. Hence, the GP description 
agrees well with the MB computation. These results conflict the findings of Ref. [17]. 

We revisited then the problem using the ansatz that incorporates a 'quadrupolar flow' 
term, used in Ref. [17]. The energy and stability of the system were again found not to be 
affected by the rotation of the trap. The resolution of this discrepancy lies in a sign error in 
the expectation value of the Hamiltonian of Ref. [17], which leads to a qualitatively different 
behavior of the properties of the system as a function of the rotation frequency Q. Our 
results are in agreement with findings in the literature for isotropic harmonic (see [16]) and 
isotropic anharmonic traps with slight anharmonicity [14, 15]. 

We should, finally, stress that the present variational approach to the stationary ground 
state is not an extensive study of the rotating attractive gas. Even though we can rule out 
the stability-enhancement of the stationary ground state and the vortex nucleation in the 
attractive gas rotating with a frequency Q smaller than the resonance frequency Q r , there is 
more physics beyond that. For instance, the stability of low- lying excited states with nonzero 
total angular momentum is expected to depend on fi. 3 We have shown elsewhere [8] that 
ground states with L > are generally fragmented states and thus more stable against 
collapse. A ground state with L = will not couple to the external rotation and hence the 

3 Indeed, we have some numerical indication for such a dependence in a MB treatment of the problem 
involving ground states of L > 0. 
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ground-state symmetry will not change as f2 increases. On the other hand, a ground state 
initially with nonzero L will be affected by the rotation. The critical parameter A c in that 
case could increase as a function of Q, before the latter reaches the extreme value Q r , and this 
will further stabilize the rotating L ^ state against collapse. It is still to be investigated 
whether and for which parameters' values crossings of energy levels and symmetry changing 
of the (noncollapsed) rotating ground-state might occur. Last, based on the found absence 
of symmetry breaking of the ground state in the examined region Q < Q r and the divergence 
of the energy and angular momentum for Q > Q r , we may speculate that, in the rotating 
attractive gas a vortex ground state - if at all can exist - may only appear as a giant-vortex 
(i.e., a single vortex at the center of the trap whose radius and vorticity are increasing 
functions of time) and this only for rotation with frequency Q > Q r . A time-dependent 
many-body treatment, for instance using the multiconfigurational time-dependent Hartree 
for bosons (MCTDHB) method [27] that has described successfully many-body dynamics of 
attractive BECs [28] , should shed light on this interesting problem and uncover the response 
mechanisms of attractive gases to rotations. 
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Appendix: The Gross-Pitaevskii energy functional with the 'Quadrupolar Flow' 
Ansatz 



We re-derive and discuss the expression for the energy functional of Ref. [17]. The GP 
energy functional in the co-rotating frame reads: 



E 



I 



h 2 



■|V</>| 2 + V(r 



\ A + ihn(<p*x^ + <p y ^-) 



dr, 



(A.l) 



2m 1 ' T 1 ' 2 1 T 1 yr dy dx 

where Ao measures the strength of the interaction, m is the mass of the particle and Q is 
the frequency of the rotation around the z-axis. In Ref. [17] there is an algebraic error in 
the above expression (Eq. (4) of Ref. [17]). There, the sign of the last term of the integrand 
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is a minus instead of a plus. This sign error remains further in the calculations of Ref. [17] 
and is seen in Eq. (8) and (10) therein. Indeed, in the last term of Eq. (10) of Ref. [17] the 
'-2' term has to be omitted and so the corrected expression would read: 

, 2 



1 

£G= 4 



. Ix ly iz 



V^lxlylz 4 ll + ly 

(A.2) 

(for the a > branch), where 7^ = l x>y / ^/h/mco and k = |Ao|iV/47r. The same correction 
is required for Eq. (8) of Ref. [17] as well. The presence of this extra term gives rise to an 
artificial dependence of the critical value of A c on the rotation frequency Q, qualitatively 
different from the correct one. Indeed, a first order expansion of the (correct) energy, 
Eq. (A.2), around 7^ = 7^, i.e., for small deformations, will result in an expression of 
the energy that does not depend on the frequency Q. According to this, for zero or small 
ellipticity e of the trapping potential, the resulting shape of the orbital (f> is symmetric around 
the z-axis, i.e., 7^ = j y , and the energy of the system, as well as the critical interaction 
strength, practically do not depend on the frequency Q. On the other hand, the (incorrect) 
energy e G as it is calculated in Ref. [17] strongly depends on Q. Furthermore, it can be 
easily seen that the 'quadrupolar flow' ansatz of either Eq. (6) or (7) gives, for small e, an 
expectation value of almost zero angular momentum (L z ) = | ^ + ^ raNVt e =° 0, and hence 
cannot describe any state with nonzero angular momentum that can in principal increase 
the stability of the system. For zero or small e the energy and the critical parameter A c 
cannot change as a function of Q and this reflects the cylindrical symmetry of the ansatze 
used, since l x ~ l y if e ~ 0. 
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FIG. 1: (Color online). Many-body calculations for the critical parameter A c as a function of the 
frequency of the external rotation O, for the cases of isotropic [e = 0, £ = 1; lower (red) line] and 
slightly anisotropic [e = 0.1, £ = 1; upper (blue) line] confining traps. The critical interaction A c 
remains practically unaffected (note the scale!) for the whole region of € [0, O r = yl — e). The 
number of particles is N = 12. See text for more details. All quantities are dimensionless. 







C=o 


C=i 


C=5 


e=0 


9.547(7) 


8.425(9) 


5.522(8) 


£=0.1 


9.554(9) 


8.429(0) 


5.523(2) 



TABLE I: Critical parameter A c for different values of the anisotropy e and the z-deformation £ 
of the trapping potential calculated in the GP theory, using the 'quadrupolar flow' ansatze, see 
Eqs. (6)-(7) and text below it. The parameter A c in the radially symmetric (e = 0) cases does not 
depend on the frequency Q, of the rotation, while the change in A c as f2 varies from (for which 
A c values are collected in the table) to Q r = yl — i is negligible (less than 0.2%) when the trap is 
slightly anisotropic (e = 0.1). All quantities are dimensionless. 
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